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ABSTRACT 

We report on the imaging analysis of ~200 ks sub-arcsecond resolution Chan- 
dra ACIS-S observations of the nearby Seyfert 1 galaxy NGC 4151. Bright, 
structured soft X-ray emission is observed to extend from 30 pc to 1.3 kpc in the 
south-west from the nucleus, much farther than seen in earlier X-ray studies. The 
terminus of the north-eastern X-ray emission is spatially coincident with a CO 
gas lane, where the outflow likely encounters dense gas in the host galactic disk. 
X-ray emission is also detected outside the boundaries of the ionization cone, 
which indicates that the gas there is not completely shielded from the nuclear 
continuum, as would be the case for a molecular torus collimating the bicone. In 
the central r < 200 pc region, the subpixel processing of the ACIS data recovers 
the morphological details on scales of <30 pc (<0.5") first discovered in Chandra 
HRC images. 

The X-ray emission is more absorbed towards the boundaries of the ionization 
cone, as well as perpendicular to the bicone along the direction of a putative torus 
in NGC 4151. The innermost region where X-ray emission shows the highest 
hardness ratio, is spatially coincident with the near-infrared resolved H2 emission 
and dusty spirals we find in an HST V — H color image. The agreement between 
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the observed H2 line flux and the value predicted from X-ray-irradiated molecular 
cloud models supports photo-excitation by X-rays from the active nucleus as the 
origin of the H2 line, although contribution from UV fluorescence or collisional 
excitation cannot be fully ruled out with current data. The discrepancy between 
the mass of cold molecular gas inferred from recent CO and near-infrared H2 
observations may be explained by the anomalous CO abundance in this X-ray 
dominated region. The t otal mass derived from the X -ray observation agrees 



with the measurement in IStorchi-Bergmann et al.l (j2009[ ) . 



Subject headings: X-rays: galaxies — galaxies: Seyfert — galaxies: ISM — ISM: 
jets and outflows — galaxies: individual (NGC 4151) 



Introduction 



NGC 4151 (classified as (R )SAB(rs)ab; Ide Vaucouleurs et al.lll99lh . one of the nearest 
{d ~ 13.3 Mpc, 1" = 65 pc; Mundell et al. 1999), and apparently brightest, active galaxy 
(Seyfert 1.5, Osterbrock & Koski 1976), has been extensively studied in almost all wave- 
bands (see Ulrich 2000 for a review). Because of its proximity, NGC 4151 offers a unique 
opportunity to examine the fueling of the central supermassive black hole (SMBH) and its 
impact on the host galaxy. 

NGC 4151 is close to face-on (inclination i ^ 21°) with a major axis position angle 
(P.A.) ^ 22° (Pedlar et al. 1992; Mundell et al. 1999). A biconical extended narrow line 
region (ENLR) is elongated along the northeast (NE) and southwest (SW) of the active 
galactic nucleus (AGN), at P.A. ~65° and 230°, respectively (Perez et al. 1989; Perez- 
Fournon & Wilson 1990; Penston et al. 1990; Evans et al. 1993). Clumpy ionized gas is 
seen in Hubble Space Telescope (HST) narrow-band images (e.g., [OIII] A5007A, Boksenberg 
et al. 1995; Winge et al. 1999) outflowing along the bicone (Schulz 1990; Crenshaw et al. 
2000; Kaiser et al. 2000; Das et al. 2005), and extending to ~1 kpc away from the AGN. The 
morphology, kinematics, and line ratio diagnostics from optical studies of NGC 4151's ENLR 
(e.g., Robinson et al. 1994; Hutchings et al. 1998; Das et al. 2005) strongly support that 
AGN radiation in the optical bicone intersects with the host galactic dis k, and photoionizes 



ambient galactic gas that is part i cipating in normal g alactic rotation f lPedlar et al.l Il992 
Evans et al.lll993l : lAsif et~aDll998[ iMundell et aDll999t . 



Near- infrared (NIR) emission line mapping of the narrow line region (NLR) of NGC 4151 
(e.g., Pa/3: IStorchi-Bergmann et al.ll2009l . l2010l ) shows that the ionized gas closely follows the 
optical bicone, and extends outside the bicone (also found in the HST study by Kraemer et al. 
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2008). The molecular hydrogen (H2) emission distribu tion is perpendicular to the bicone axis 
( iFernandez et al.l Il999l : IStorchi-Bergmann et al.l 120091). following approximately the P. A. o f 
a large scale (3.2^ x 2.50 g as-ric h stellar bar f^l30°: IPedlar et al.lE99^ : iMundell et al.lll999[ ). 
Storchi-Bergmann et al.l (120101 ) propose that H2 emission is potentially excited by X-rays 
from the AGN or shocks related to t he accretion flow al ong the large scale bar. A spatially 
resolved study of the ^^CO emission (IDumas et al.ll2010l ) finds two prominent CO gas lanes 
at 1 kpc from the nu cleus, but no cold rn olecular gas in the central 300 pc where the NIR H2 
line emission peaks. iDumas et al.l (120101 ) further suggest that the H2 line emission is likely 
photo-excited by the AGN. 

The ENLR of NGC 4151 has been resolved in the soft X-rays with Einstein and ROSAT 
observations (Elvis, Briel, & Henry 1983; Morse et al. 1995). Early Chandra images show 
that the extended soft X-ray emission is generally correlated with the optical forbidden line 
emission (e.g., Ogle et al. 2000; Yang et al. 2001), whereas the nucleus remains unresolved in 
the hard X-ray band. Our recent deep Chandra observations of NGC 4151 (PI: Fabbiano) on 
the circum-nuclear region are now providing detailed information on several spatial scales. 
The Chandra HRC image (Wang et al. 2009) resolves the nuclear X-ray emission on spatial 
scales of ~30 pc, finding evidence for interactions between the radio outflow and NLR cloud^. 
Our deep ACIS data have led to the discovery of a large scale soft diffuse X-ray emission 
(Wang et al. 2010a) extending out to ~2 kpc from the active nucleus and filling in the 
cavity in the HI 21 cm emission distribution (lAsif et al.lll998l ; iMundell fc Shondll999l ). This 
extended X-ray emission could be due to either hot gas in the galactic disk heated by the 
nuclear outfiow or gas photoionized by an outburst of the AGN (Komossa 2001; Wang et al. 
2010a). 

This is the first of a series of papers where we report on the results from deep Chandra 
ACIS observations of the circum-nuclear region in NGC 4151, focusing on the imaging anal- 
ysis and relating the X-ray morphology to the cold interstellar medium (ISM) in the host 
galactic disk. This paper is organized as follows. In § 2 we describe the observations and data 
reduction; in § 3 we describe the imaging analysis and present the X-ray morphology. In § 4 
we discuss the results, in particular the role of the X-ray emission in exciting the molecular 
material near the AGN. We summarize our findings in § 5. In a forthcoming paper (Wang et 
al. in preparation; Paper II), we will present the X-ray emission line strength maps, detailed 
comparison of these maps with the ionized gas distribution using Ha and [OIII] images, and 
spatially resolved spectral fitting with self-consistent photoionization models. 



^For a comprehensive view of the NGC 4151 nuclear structure, see Mundell et al. (2003). 
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2. Observations and Reduction 
2.1. Chandra data 

NGC 4151 was observed with Chandra for a total of ~200 ks, using the spectroscopic 
array of the Advanced CCD Imaging Spectrometer (ACIS-S; Garmire et al. 2003) in 1/8 
sub-array mode during March 27-29, 2008. The nucleus of NGC 4151 was placed near the 
aimpoint on the backside-illuminated S3 chip, which has a high sensitivity to soft X-rays. 
The data were analyzed following the standard procedures using (Version 4.2) with 

CALDfllfl (Version 4.2.1) provided by the ChandraX.-r&y Center (CXC). The lightcurve from 
regions free of bright sources on the S3 chip was used to remove times of high background 
count rates. The cleaned data have total good exposure times of 180 ks (116 ks and 63 ks, 
in ObsID 9217 and ObsID 9218, respectively). Table 1 summarizes the observations used in 
this paper. 

On the 1/8 S3 chip, a total of 24 point sources were detected using a significance 
threshold 10~^ with the wavelet-based algorithm wavdetect (Freeman 2002), part of the 
CIAO tools. The search was performed in wavelet scales from 1 to 16 pixel (0.5" to 8") 
in steps of \/2. These sources are removed from the images when we perform the spectral 
extraction of the diffuse emission. In addition, matching the positions of the X-ray detections 
to the SDSS Data Release 7 catalog indicates an absolute astrometric accuracy of 0.3". 

The two ACIS observations of NGC 4151 were taken almost uninterrupted, with identical 
roll angles. The events were reprojected to the reference frame of ObsID 9217, and merged 
to create a single image. The pipeline randomization was removed with the CIAO tool 
acis -process -events. Subpixel event repositioning and binning techniques (Tsunemi et al. 
2001; Mori et al. 2001; Kastner et al. 2002; Li et al. 2003, 2004) were applied to improve 
the effective resolution of the ACIS images to better than 0.4" (Li et al. 2004). The images 
with subpixel resolution presented in this paper used the energy-independent shift for the 
corner split and 2-pixel split events ("static" method in Li et al. 2004). 

Spectra and instrument responses were generated using CIAO Version 4.2. The source 
spectra were analyzed with XSPEC Version 12.5 (Arnaud 1996), taking background spectra 
from a blank region from the same CCD node outside of the galaxy. Spectra were grouped 
to have a minimum of 20 counts per energy bin to allow for fitting. When extracting 
spectra, the BACKSCAL keyword (equal to the fractional area of the chip occupied by 



^ ^http : //cxc . harvard . edu/ciao/ 
" http : //cxc . harvard . edu/ caldb/ 
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the extraction region) was computed for the extended regions using CIAO tool s^^eceztracfl. 
We restricted our modehng to photon energies between 0.3 keV and 7 keV, to avoid the 
instrument cahbration uncertainty at lower energies, and the higher background at high 
energies. The 90% confidence interval for a single interesting parameter is reported for all 
spectral fitting results. 



2.2. Archival Optical/NIR imaging 



To examine extinction with a color map of the circum- nuclear region of NGC 4151, we 
have also retrieved from the Hubble Legacy Archiv^the HST ACS-HRC/F550M (\^-band) 
and NICMOS-NIC2/F160W (if-band) images. These images were acquired as part of the 
Cycle 12 Program 9851 (PI: Peterson) and C ycle 7 Program 7215 (PI: Thompson), and the 
details of these observations were published in lOnken et al.l (120071 ) and [Thompson fc Corbin 
( 1l999l ). respectively. Pipeline processed images are further cleaned to remove cosmic rays 
using the Laplacian Cosmic Ray Identification algorithm (van Dokkum 2001). Model PSFs 
were generated with the TinyTim tool (Krist 1995) for the specific instrument and filter, and 
subtracted to remove the unresolved nuclear component. The resolution of the ACS l^-band 
image has been degraded to match that of the if-band, by convolving with the broader 
NICMOS PSF {FWHM = O.IQ"). 



3. Imaging Analysis and Results 

3.1. The "Raw" ACIS Image 

Chandra ACIS images of the central ~ 30" x 30" (~ 2 x 2 kpc) region of NGC 4151 
were extracted in the 0.3-1 keV and 1-7 keV energy bands based on its X-ray spectrum. 
These images, shown in Figure [H are in good agreement with the ACIS-S image obtained by 
Yang et al. (2001) and the HETG zero-order image shown by Ogle et al. (2000). However, 
our deeper images are ~7 and 20 times more sensitive than the previously published ACIS 
images of this circum-nuclear region, respectively. Figure [1] clearly shows the bright nucleus 
dominating at energies >1 keV and soft resolved extended regions (the streaks on both sides 
of the nucleus are an artifact of the CCD readout process). 



http : //cxc . harvard . edu/ ciao/ threads/specextract/ 
http : //hla . stsci . edu 
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The SW side (along P.A.~233°) of the 0.3-1 keV emission is clumpy and clearly seen to 
r = 14", with fainter diffuse emission extending as far as 20" from the nucleus. A relatively 
bright "blob" to the SW, noted in Yang et al. (2001), is located at r = 6" (380 pc) from the 
nucleus at a P. A. of ~227° (Figure [H indicated as "A" with an arrow). Our image also shows 
at least one more such clump of X-ray enhancement further out at 10" from the nucleus ("B" 
in Figured]). On the NE side, the soft emission is brightest at r < 5" and extends to 11" 
from the nucleus along P.A.~75°. The soft X-ray emission along the NW-SE directions is 
fainter than the NE-SW emission by a factor of 2 and is less extended (3"; ~200 pc). 

3.2. The Nucleus: Photon pile-up and the Nuclear Point Spread Function 

(PSF) 

The C/ian(ira-measured position of the Seyfert X-ray nucleus is ax(2000) = 12^^10^32. 59*^, 
(5x(2000) = +39°24'21.2". Given the uncertainty in the absolute astrometry, this position is 
identical to the optical (a = 12'^10"'32."6, 6 = +39°24'21", Clements 1981) and radio nuclear 
positions {a = 12'^10™32.*6, 5 = +39°24'21", Ulvestad et al. 2005). 

The nucleus of NGC 4151 is a well known strong X-ray source. During the exposure of 
a single ACTS frame, two or more photons from such a high count rate X-ray source may 
be recorded as a single event. This is known as the "pile- up effect"|^, which causes the loss 
of information from the original events, resulting in a hardened source spectrum and grade 
migration (Davis 2001). Aiming to study the extended X-ray emission as close to the nucleus 
as possible, our observations were performed using the ACTS "1/8 sub-array" mode (0.6 s 
frame time) to reduce nuclear pile-up. However, pile-up is still a concern as Yang et al. 
(2001) showed that, although their observation was taken when NGC 4151 was experiencing 
a low flux state, pile-up still significantly affects the inner 1.3"-radius region in the ACIS 
image in the same mode. Before we proceed to further imaging analysis, we must evaluate 
the effects of pile-up on the nuclear PSF, to exclude that some of the observed features may 
be due to pile-up effects. 

We first compared the radial profiles of the observed X-ray emission to that expected 
from a point source. To simulate the nuclear PSF, we need to know the intrinsic spectrum 
of the heavily piled-up nucleus. The events extracted from the ACIS readout streak do not 
suffer pile-up as the incoming photons during each frame readout were still recorded by the 
ACIS CCD with an exposure equivalent to 40 fis per frame. We used these readout events to 



^For more information see Chandra ABC Guide to Pile Up, available at 
http : / / cxc . harvard . edu/ciao/ download/ doc/pileup_abc . ps 
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model the spectrum of the nucleus (Wang et al. 2010b) and obtained an absorption corrected 
luminosity of L2_iokcV = 1-3 x 10"^^ erg s~^. 



Carter et al. 2003) which was run with the same observation configurations as our ACIS 
observations (Table 1). The ChaRT output was then fed into the MARX simulatoj^ to 
project rays onto the detector and create an event file, with responses applied. We also used 
MARX to simulate a piled event file using the ~35% pile-up fraction found by Wang et 
al. (2010b), in order to compare with the observed radial profiles. The effectiveness of this 
procedure was demonstrated by Russell et al. (2010). 

Figure [2^ compares the observed radial profiles of the nuclear region of NGC 4151 in 
the 1-7 keV band with the simulated PSF of a point-like nucleus (with and without pile- up). 
The count rate in the brightest inner 2 arcsec region is significantly lower than the expected 
values without pile-up, and consistent with the PSF including MARX simulated pile-up. The 
soft (0.3-1 keV) X-ray emission close to the nucleus is clearly resolved, both along the NE - 
SW and the NW - SE direction. In both directions, the soft extended emission dominates 
over the nuclear PSF wings at radii > 2" from the nucleus (Figure [2)d and[2t). Beyond a 
radius of ~3" from the nucleus, the piled MARX simulations produced essentially identical 
results to the simulations without pile-up, so at these larger radii pile-up is not an issue. 

As pile-up combines events of different energies and causes grade migration, another 
way to estimate the importance of pile-up is using the ratio of events with "bad" ASCA 
grades (1, 5, 7) to those with "good" ASCA grades (0, 2, 3, 4, 6) where a bad/good grade 
ratio of < 0.1 indicates that pile-up is approximately < 10% (Russell et al. 2010). The radial 
distribution of the bad/good grade ratio is shown in Figure [2]i for both the 0.3-1 keV and 
the 1-7 keV bands, which indicates that, beyond r ~ 2" from the nucleus pile- up is no longer 
significant in the 0.3-1 keV band (and beyond r ~ 3" for the 1-7 keV band). We conclude 
that for our analysis of the extended soft X-ray emission, the effect of pile-up is mild beyond 
a radius of 2" from the nucleus, and effectively absent at r > 3". 

For imaging analysis, the major effect of pile-up on the ACIS on-axis PSF is that, the 
core (where most of the photons hit) gets suppressed, because some piled events are rejected 
as bad events, creating a circular "caldera" (e.g., Yang et al. 2001). At the pile-up fraction 



^High Resolution Mirror Assembly (van Speybroeck et al. 1997); see 
http : //cxc . harvard . edu/ cal/Hrma/u sers_guide/, for details. 

^ http : //cxc . harvard . edu/ chart/ 

' http : //space .mit . edu/ cxc/marx/ 



This model was used as the input spectrum to the Chandra HRMA|3 ray tracer 
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of ~35% encountered in this observation, there is ~40% decrease in brightness in the core of 
the PSF, but httle distortion of its circular shape (Gaetz et ah, private communication; see 
the Chandra Proposers' Observatory Guide [POGJ^ Section 6.15.). Thus, ahhough we will 
discuss the images of extended features in the central 2" region that are affected by pile-up 
with caution, deviations from circular shape must be related to real X-ray emission features. 



3.3. ACIS Images with Subpixel Binning 

Chandra HRMA's PSF (half-power diameter HPD=0.6", E < 1.5 keV; see the Chandra 
POG) is undersampled by the ACIS CCD pixels (0.492" x 0.492"). To take advantage of the 
telescope dithering^ to improve the sampling, a fine pixel size (0.0625", ~l/8 of the native 
ACIS pixel size) was used when extracting the images. This subpixel binning approach is 
frequently adopted in imaging studies of X-ray jets pushing for the highest spatial resolution 
(e.g., Harris et al. 2004; Siemiginowska et al. 2007; Perlman et al. 2010). In this section we 
will also make use of the 50 ks Chandra observation of NGC 4151 with the High Resolution 
Camera (HRC; Murray et al. 1997) presented in Wang et al. (2009). 

Figure |3] compares the 0.3-1 keV subpixel ACIS and HRC (Wang et al. 2009) images both 
smoothed with a FWHM = 0.4" gaussian kernel. Although the raw soft band ACIS image 
(Figure [2^) already shows the presence of extended emission along the NE-SW direction 
and hints of split emission in both NE and SW cones, the subpixel processed ACIS image 
of Figure [3^ reveals more morphological details, especially the presence of several X-ray 
blobs (marked as 'A" , "B" , "C" , and "D" in Figure |3]) and the clear bifurcation of both 
the NE and SW cones. The appearance of the NE-SW soft X-ray emission is an elongated 
"X" -shape with an acute apex angle, resembling edge-brightened cones (c.f. Figured^). We 
have measured the P.A. of the brighter bifurcated SW X-ray emission: the edges are located 
at 227° ± 8° and 245° ± 5°, with a full opening angle of ~ 20°. The HRC image shows 
the identical extended morphology for the circum- nuclear region (Figure [3]d), validating that 
subpixel repositioning does not introduce unwanted artifacts. 

Figures and b zoom in on the inner ~ 6" x 6" (390 pcx390 pc) nuclear region, with 
and without subpixel resolution. The X-ray emission on scales of 200 pc (~ 3") is clearly 
extended along the NE-SW direction. This extent is in good agreement with the overall 
morphology of the [OIII] emission (e.g. Kaiser et al. 2000; Wang et al. 2009). The ACIS 



^'^ Available at http : //cxc . harvard . edu/proposer/PDG/ 
"see POG Version 12 § 5.8.2, page 75 (2009). 



- 9 - 



images with subpixel binning (e.g., Figures [3^ andHb) reveal curvy extensions 2" away from 
the nucleus (labelled as "C" and "D"), which are not discernable in the raw soft ACIS image 
(Figure Sh-). 

The subpixel ACIS image also shows a linear feature at P.A.=235° (labelled as "E" in 
Figure Hb) and discrete knots in the inner region (on scales <1" from the nucleus). The soft 
X-ray morphology is similar to that of the NLR clouds seen in HST optical images of NGC 
4151 (e.g., Evans et al. 1993, Winge et al. 1999; see Wang et al. 2009 for a Chandra/RRC 
and HST/WFPC2 image comparison). Corresponding outer clouds were previously detected 
in deep HST images (Hutchings et al. 1998; Kaiser et al. 2000). Comparison of these features 
to the [OIII] emission line clouds and radio outflows (e.g.. Pedlar et al. 1992; Mundell et al. 
2003), and emission modeling will be presented in Paper II. 

We emphasize that although the inner 2" region is affected by pile-up (§ 13. 2p in the 0.3-1 
keV band, the ACIS features seen in the soft X-ray emission are not observation speciflc or 
pile-up artifacts, for the following reasons: 



For the on-axis nucleus, the PSF with pile-up simulated with MARX is circular, thus 
the elongated features (4" across) seen in the soft band cannot be due to pile-up related 
PSF distortion; 

Using only the single pixel events (grade 0), which give the best spatial resolution and 
are likely to be unaffected by pileup (Ballet 1999; Davis 2001), the extended features 
are clearly present (Figure Hb); 

Low pile-up ACIS images from Chandra archival observations of NGC 4151 (ObsID 
372) show structures of similar shape and spatial scale as the features seen in our 
images (Figures SJi and e). Because of the short frame time (0.1 s and 0.4 s) and 
a factor of two lower nuclear flu x, the estirnated pile-up fraction is 6% a nd 11% in 



these observations, respectively (lYang et al.ll200ll : iGonzalez-Martfn II2008I ). We also 



examined the HETG zero-order images of NGC 4151. Although some suffer signiflcant 
pile-up due to longer frame time and a brightened nucleus (ObsID 7830), extended 
emission along NE-SW with similar substructures as we detect is seen in ObsIDs 335 
and 7829; 

The PSF- deconvolved HRC image (Figure 2 in Wang et al. 2009), which is not affected 
by pile-up, shows the same extended morphology and small-scale enhancements within 
2" from the nucleus that are not seen in ACIS without subpixel binning. 



We applied two image reconstruction techniques to the ACIS data, the Richardson-Lucy 
(R-L) algorithm (Richardson 1972; Lucy 1974) and the Expectation through Markov Chain 
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Monte Carlo (EMC2) method (Esch et al. 2004; Karovska et al. 2007). Both show similar 
extended morphology, although the R-L deconvolution produces more grainy restoration of 
extended features. Figure O compares the EMC2 PSF deconvolved ACIS image with the 
similarly deconvolved HRC image from Wang et al. (2009). The resemblance of the NE-SW 
elongation and substructures (such as the curvy feature "C" and NE clump "D" marked in 
Figure |3^) between the subpixel ACIS images and the deconvolved HRC strengthens our 
confidence that the resolved features are real. However, we noticed that the X-ray clump 
1" NE of the nucleus in the ACIS image is displaced in the HRC image (~0.3" closer to 
the nucleus). The morphology of the enhancement located 1" SW to the nucleus is also 
different in two reconstructed images: it is mostly a linear north-south elongation in the 
ACIS image, whereas in the HRC image it shows a tail towards SW). Two factors may cause 
these discrepancies: (1) The region within l"radius of the nucleus in the ACIS image is 
affected by pile-up, but the accuracy of the piled PSF on this scale has not been calibrated; 
yet this inaccuracy may have an impact on the reconstructed image; (2) HRC is sensitive to 
photons in the energy range of 0.1-10 keV, therefore may reveal real structures that are not 
seen in the 0.3-1 keV ACIS image. 



3.4. Hardness Ratio Image of the Central 3 kpc-diameter Region 

To find regions of harder spectral index or higher X-ray obscuration, we extracted 
Chandra images in the 0.3-1.0 keV and 1.0-2.0 keV bands and derived the counts ratio 
Ci-2kev over Co.s-ikev- These bands are sensitive to hydrogen column densities of the order 
of 10^^ cm~^, and the relative changes in the ratio map have been demonstrated to reflect 
variations in the X-ray absorption (Bianchi et al. 2007). The resulting hard ness ratio map 



of th e central 50^^ (3.2 kpc) diame ter region is shown in Figure E^, with Ha (IKnapen et al. 



20041 ) and CO (jPumas et al.ll2010l ) contours outlining the ionized gas and the cold molecular 
gas lanes, respectively. Figure [6)d shows the hardness ratio map of the inner 6" (400 pc) 
diameter region. 

One obvious feature in Figure [6^ is the harder (white) region extended in the NW-SE 
direction (P. A. ~150°). This is approximately the P. A. of the minor axis of the galaxy and 
the weak bar (Mundell & Shone 1999), and perhaps also t he dir ection of a putative nuclear 



torus, mapped in H2 emission by IStorchi-Bergmann et al.l (120091 ) (contours in Figure Eb; see 



also Fernandez et al. 1999). Thus higher obscuration is plausible in this region. Coincident 



Asif et al. 


1998; 


Ohtani et al. 


2001 


), noted by 


Terlevich et al. 



their V — R image as "the red bar" (P.A.=150°; 1" in length). 
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To examine the extinction of the nuclear region of NGC 4151, we created the V — H color 
map using HST images (§ 12. 2p . covering the central 10" x 10" of NGC 4151 (Figure [7]). Dark 
spiral features of higher extinction are clearly revealed towards the nucleus perpendicular to 
the bicone directio n, thanks to the high spatial resolution afforded by HST. Adopting the 



extinction curve in ICardelli et al.l (119891 ) . the observed E{V — H) suggests a range between 
0.5 and 1 for E[B — V), which is equivalent to 1.5 < Ay < 3 (mag). 

Another notable feature in Figure is that the edges of the ENLR show harder emis- 
sion than the inner region. To obtain quantitative constraints on the hardness ratios, we 
calculated the uncertainties using the Bayesian estimation of Hardness Ratios (BEHR; Park 
et al. 2006) tool provided by the California-Harvard Astro-Statistics Collaboration!^, where 
the hardness ratio, HR, is defined in BEHR as (Ci_2kcv — C'o.3-ikev)/(C'i-2kcV + C'o.s-ikev)- 
For the inner part of the cone (named Rl) HR = — 0.63!1]q'q25 edges of the optical 

outflow (R2) HR = —0.27 ± 0.03, and for the highest hardness ratio region perpendicular 
to the cone-axis (R3) HR = —0.04 ± 0.02. Thus the hardness ratios differ significantly for 
each of these regions. This is consistent with the hollow bi-cone geometry proposed in Das 
et al. (2005) (see also Evans et al. 1993; Robinson et al. 1994; Hutchings et al. 1998; 
Storchi-Bergmann et al. 2010), where clumpy outflowing material is mostly located along 
the boundaries of the bi-cone. One side of the cone is expected to intersect an area of the 
host galactic disk, and the projected appearance will be an illuminated narrow sector (e.g., 
region Rl) with higher absorption along the edges, which is indeed observed here (region 
R2). 

To further investigate these absorption column density variations, spectra from the three 
regions above were extracted using the hardness ratio map, as outlined in Figure [HI Because 
the PSF scattered emission is not negligible for the bright nucleus, we fit the spectra with 
a simple model consisting of an absorbed soft thermal component and a nuclear component 
(Wang et al. 2010b), allowing the normalization of the nuclear component to vary to account 
for contribution from the PSF wings. Figure [9^-c show the spectra and the fits; the fitting 
results are summarized in Table [2l These spectra confirm the wide range of X-ray absorption 
suggested by the hardness ratio map. The edges of the cone (R2) indeed show a higher 
absorbing column [Nh = 7± 5 x 10^^ cm~^) than the inner cone region (Rl), for which only 
a Galactic column {Nfj = 2x 10^*^ cm~^; Murphy et al. 1996) is required. The highest column 
density, Nh = 6.5 ± 1.5 x 10^^ cm~^, is found across the nucleus (R3) and is approximately 
equivalent to E(B - V) = 10 adopting N{HI + H2)/E{B - V) = 5.8 x 10^^ cm'^ mag"^ 



(IBohlin et al.l 119781 ). This value is significantly higher than the average E{B ~ V) ^ 0.5 
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found in the NL R (IStorchi-Bergmann et al.ll2009l ) and the E{B — V) = 0.5 — 1.0 found in 
the V — H map. iMundell et aLl (119951 ) find an HI absorption column density of 3.9 x 10^^ 
cm~^ towards the radio nucleus; the factor ~15 higher X-ray Nh suggests that most of the 
intervening hydrogen is in the form of H2 or high abundances ~ 15Z/Zq. 



4. Discussion 

4.1. Extended X-ray Morphology 

This is the first time such a high sensitivity high resolution X-ray spectral imaging 
of the NGC 4151 nuclear region has been obtained. The HRC image provides no spectral 
information, and previous ACIS images (Ogle et al. 2000; Yang et al. 2001) either lack 
sensitivity to convincingly identify the X-ray features seen in subpixel resolution, or suffer 
from heavy pile-up in the inner 3". 

The morphological features of the X-ray emission of NGC 4151 are summarized as 
follows: (1) By pushing the spatial resolution of the ACIS image with subpixel resolution 
techniques. X-ray enhancements in the innermost nuclear region (r < 1") seen in the PSF- 
deconvolved HRC image are recovered; (2) The extended soft X-ray emission shows the same 
biconical morphology as the optical ENLR out to a radius (20"; ~1.3 kpc) larger than seen 
in earlier X-ray studies (c.f. 14"; ~900 pc, Yang et al. 2001); (3) There is also firm detection 
of extended X-ray emission at angles beyond the boundaries of the optical ENLR ionization 
cone; (4) The observed X-ray emission is more absorbed towards the edge of the cone and 
perpendicular to the cone axis at the apex of the cone. 

The two-sided, bifurcated X-ray emission morphology is c onsistent with the model de- 



scribed in Das et al. (2005) and lStorchi-Bergmann et al.l ( I2OIOI ). in which clouds are outflow 



ing along the walls of a hollow bicone inclined 45° to the line of sight. One side of the cone 
is expected to exit the plane at 21°, interse cting an elongated area of the host galactic disk. 



The projected appearance (c.f. Figure 9 in IStorchi-Bergmann et al.ll2010l ) closely resembles 
the observed morphology here. 

The extended X-ray emission seen perpendicular to the ENLR, along the NW-SE direc- 
tion, was first suggested to be "unresolved or marginally resolved" in Yang et al. (2001). The 
definitive resolved emission at angles beyond the edges of the bicone (Figure ^) indicates 
that the gas there is not completely shielded from the continuum, as would be the case for 
a molecular torus collimating the bicone. This agrees well with the detection of significant 
optical emission outside the emission-line bicone in the HST [OIII] and [OH] images, where 
Kraemer et al. (2008) were able to model the emission lines with a weaker ionizing con- 
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tinuum, filtered by an io nized absor ber partially covering the nuclear continuum, i.e. 
accretion disk wind 



the 



le.g., 



Elvisll2000h . 



The NE soft X-ray emission reaches as far as 11" (700 pc) from the nucleus, then the 
surface brightness decreases rapidly and significantly, by a factor of 5. This location lies 
near to both the northern CO g as lane and one of the dust arcs identified in the optical 
V — I color map (lAsif et al.lll998l . a region not covered by our V — H map). These features 
are gaseous compressions in ci rcumnuclear gas spiral arms due to the g ravitational influence 
of the larg e-scale stellar b ar (IMundell et al.l 1 19991 : iDumas et al.l |2010| ) . Since the biconical 
geometry (jPas et al.ll2005l ) implies that both cones are intersecting the host galactic disk, 
the terminus of the NE X-ray emission is likely to occur where the NE outflow encounters 
dense ISM in the galactic disk. We will discuss this boundary further in Paper II. 



4.2. X-ray Absorption Features and Cold Material for AGN Feeding 



The feat ures seen in the hardness ratio map (Figure 



cone mode l s (jPas et al.ll2005l : IStorchi-Bergmann et al.ll2010l ). In the geometry illustrated by 



6^) also agree well with the bi- 



Das et al.l (120051 ) . the more absorbed X-ray emission seen along the bi-cone edges (R2 in 
Figure [S]) can be explained by higher absorption from the clumpy medium outflowing along 
the wall of the SW cone, as it emerges from the galactic disk. Adopting the conical geometry 
in Das et al. (2005) with an inner full opening angle of 30° and outer full opening angle 
of 66°, and assuming that the Nh is due to the extra optical path along the cone wall, we 
estimate an average uh ~ 10 cm~'^ in the cone wall (c.f. nu = 220 cm~^ for optical clouds; 
Penston et al. 1990). The region perpendicular to the bi-cone axis (R3 in Figure [HD may 
represent the inward extension of the CO gas lanes. 

The innermost r < 1.5" region (P.A.~150°) of high X-ray har dness ratio (Figure Eb) is 
spatially coincident with the H2 A2.1218 /xm emission (grey contours: IStorchi-Bergmann et al 
2OO9I : IPernandez et al.lE999l ) and the dusty spirals in the V — H color map (Figure [7]), which 
may trace infl owing material feeding th e AGN . Next we discuss three mechanisms that were 
considered by IStorchi-Bergmann et al.l (120091 ) for the origin of H2 line emission: (1) fluores- 
cence by FUV photons; (2) thermal excitation by shocks in the accretion flow towards the 
nucleus, or (3) photo-excitation by X-rays from the AGN. 



(1) Fluorescence by FUV photons is ruled out by IStorchi-Bergmann et al.l (120091 ) based 
on the H2 2-1 S(l)/l-0 S(l) line ratio and strong evidence for thermal equilibrium {Texc = 
2155 K; IStorchi-Bergmann et al.l 120091 ) in the H2 lines. However, it is worth noting that 
the peak density of 10^-10^ cm~3 in the H2 emitting region, taken at the measured value 
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within the inner 0.5 arcsec in Storchi-Bergmann et al. (2009), is above the critical density 
of the rotational- vibrational transitions of H2 (~ 10^ cm~'^; Tielens 2005). It implies that 
thermalization by collisions effectively erases the signature of fluorescent excitation (Black 
& van Dishoeck 1987; Draine Sz Bertoldi 1996; HoUenbach & Tielens 1999), therefore the 
line ratio diagnostics used at lower densities may become uncertain to distinguish whether 
UV fluorescence or collisional excitation is dominating the H2 emission. 

(2) H2 excitation by shocks in the accretion flow is not favored by current data. Mundell 
& Shone (1999) interpreted the H2 emission as being associated with the molecular torus 
obscuring the AGN which would explain its orientation perpendicular to the radio jet. The 
inflows observed in HI observations (Mundell & Shone 1999) are only traced to the inner part 
of the oval bar (r ^ 20"; Asif et al. 1998). On the 1 kpc scale, most of the cold molecular gas 
is distributed along the curved CO lanes ~15" away from the nucleus (Dumas et al. 2010) 
with kinematics consistent with the disk and the bar (e.g., Mundell et al. 1999). No strong 
kinematic shocks are expected, which is also consistent with the gravitational potential that 
shows no important non-axisymmetric component on the relevant scales. Moreover, the 
detection of H2 in this relatively high density region {uh ~ 10^ cm~^) is evidence against the 
presence of relatively strong shocks {vg > 25 km s~^) in this region, which would dissociate 
the molecular hydrogen coUisionally. However, it remains possible that weaker shock flows 
exist and contribute to the H2 excitation. 



(3) Ph o to- exc itatio n of H-} by the AGN caxi 



Riffel et al.l ( 120081 ) and IStorchi-Bergmann et al. 



)e ev aluated more quantitatively following 
2OO9L using the line intensities calculated 



for X-ray-irradiated molecular gas by lMaloney et al.l (119961 ) and our measured parameters. A 
key p arameter is the effect ive ionization parameter, ^eff = 1-26 x / {n^N 22) (Equation 

12 in iMaloney et al.lll996l ). where is the incident X-ray flux, is the cloud density in 
units of 10^ cm~^ and N22 is the absorption column in units of 10^^ cm~^. For F^ = 7.6 
erg cm^^ s^^ (-^2-iokeV = 1 x 10"^^ erg s~^; Wang et al . 2010b) at the H2 emission peak 
(~30 pc from the nucleus; IStorchi-Bergmann et al.l 120091 ) . 77,5 = 1, and A^22 = 6.5 (Table 2), 
we obtain ^efi ~ 2 x 10~^ and an emerging H2 1-0 S(l) intensity of 



sr . For the aperture (0.3 x 0.3 arcsec or 2.16 x 10 sr) used by IStorchi-Bergmann et al. 



8 X 10 ^ erg cm ^ s ^ 



(l2009l ). the expected value is 2 x 10 erg cm ^ s ^ This is a factor of two higher, but 
likely consistent with, the obs erved line flux for the same transition (1 x 10~^^ erg cm~^ s~^; 



Storchi-Bergmann et al.l 120091 ). given the uncertainties in our assumptions and the models. 
Thus, excitation by the X-rays from AGN can rea sonably explain the observed H 2 line 
emission, in agreement with the conclusion drawn in IStorchi-Bergmann et al.l (120091 ). This 
also indicates the X-ray radiation received by the H2 gas is comparable to that observed 
along our line-of-sight, which is not blocked by a molecular torus in this direction. 
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The X-rays from the AGN may also expl ain the non-detection of CO ga s in contrast to 
the 10^ — IO^Mq molecular gas estimated by IStorchi-Bergmann et al.l ( l2010l ). Dumas et al. 
(2010) showed that no CO gas was detected within the inner 300 pc of the active nucleus, 
where strong H2 line emission arises. Using a standard Galactic CO-to -Hg conversion factor 
of Xco = 2 X 10^0 cm-2 (K km s-^)'^ (Solomon & Barrett 1991), IPumas et all J2010h 
derived a 1-a upper limit of lO^M© for the cold molecular gas mass. They suggested that 
CO itself may not be a good tracer of molecular gas in X-ray dominated regions (XDR) 
close to an active nucleus, where equilibrium molecular abundances are affected (Meijerink 
& Spaans 2005, 2007). 

Model 3 in Meijerink & Spaans (2005) calculated the CO/H2 abundance ratio for = 
1.6 erg cm~^ s~^ and n = 10^'^ cm~^ (the model parameters most close to the flux and 
density seen here), which is ^ 10^^ at column densities below Nh = 10^^ cm~^. For clouds 
closer to the nucleus or immersed in diffuse X-ray emission, can be higher, and a CO/H2 
ratio of > 10~^ can only be found at column densities Nh ^ 10^^ cm~^ (model 4 in Meijerink 
& Spaans 2005). Given that the absorption column from the spectral fit is Nh ~ 6.5 x 10^^ 
cm~^, the lack of CO perhaps can be explained by a low CO abundance relative to H2 in 
the XDR. 

Assuming that most intervening material is in the form of H2, the X-ray derived Nh = 
6.5 ± 1.5 X 10^2 indicates a total H2 mass of 2 x lOj M^ in the 3" x 3" (200 pc diameter) 
nuclear region imaged by Storchi-Bergmann et al.l (2009). This value a g rees w ell with the 
lower end of total H2 mass (10^ — 10^ Mq) in IStorchi-Bergmann et al.l (120091 ) . which was 
derived using the mass of "hot" H2 (240 M©) detected in the near-IR. 

One clear prediction of the XDR scenario is that, like in shocks, the [01] 6 3/im/[CII] 
158/im ratio is higher (> 10) than generally found in PDRs (IMaloney et al.lll996l ). This can 
be tested with new observations using far-infrared (FIR) facilities such as PACS spectrometer 
(Poglitsch et al. 2010) on board the Herschel Space Observatory (Pilbratt et al. 2010). NGC 
4151 was observed as part of Herschel Key Program SHINING (PI: Sturm). They find that 
the [O I] 63/im emission is highly concentrated on the central spaxel, while the [C II] 158/im 
emission is more extended. The ratio of [O I] 63/im / [C II] 158/im in the central spaxel is 
high, ~11 (E. Sturm, private communication), consistent with XDR or shock models with 
an effective ionization parameter of log^efi ~ —3 (Figure 12 in Maloney et al. 1996). We 
note that this value is higher than log^es = —3.7 in the XDR model analysis above using 
the observed X-ray flux and column density. This discrepancy might be attributed to more 
than one mechanism at work exciting the H2 emission. Similar [OI]/[CII] ratios of ~10 have 
been reported in dense PDRs with no significant X-ray excitation, for example in regions 
near the ionizing stars of the Orion Nebula (Boreiko & Betz 1996). 
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In short, more investigations on the importance of UV fluorescence, shock excitation, 
and X-ray excitation of H2 emission in NGC 4151 are needed to further examine the caveats 
discussed above. In particular, theoretical effort on models that focus on the astrophysical 
environment around AGN will be valuable to clarify the ambiguity of the H2 excitation 
mechanism in NGC 4151 and in other systems in general. 

5. Conclusions 

We have obtained and analyzed deep Chandra observation of NGC 4151, which demon- 
strated the imaging power of ACIS with subpixel resolution to resolve the circumnuclear 
X-ray emission of a nearby Seyfert galaxy. 

The X-ray structures in the nuclear region (inner 300 pc) seen in our HRC study of 
NGC 4151 (Wang et al. 2009) are recovered in the deep ACIS image with subpixel resolution 
techniques. Moreover, the extended soft X-ray emission of NGC 4151 reaches to a distance 
of ~1.3 kpc to the SW of the nucleus, much farther than seen in earlier X-ray studies. We 
suggest that the terminus of the NE X-ray emission is likely where the NE outflow encounters 
dense ISM in the galactic disk. Extended X-ray emission is firmly detected outside the 
boundaries of the ionization cone, which indicates that the continuum is not completely 
blocked by a molecular torus. This is consistent with the findings of optical studies where 
nuclear ionizing photons were filtered by the bi-cone wind. 

The hardness ratio image and spectral fits show that the X-ray emission is more absorbed 
towards the edges of the cone, as well as perpendicular to the cone axis along the P.A. of 
the putative nuclear torus. These regions of different spectral hardness are consistent with 
the expected features in the biconical outflow models. 

X-ray emission shows the highest hardness ratio in the innermost region, spatially co- 
incident with the near-infrared resolved H2 emission and the dusty spirals revealed in the 
HST {V — H) color map. All three features possibly trace the cold material feeding the 
central SMBH. The agreement between the observed H2 line flux and the predicted value 
from X-ray-irradiated molecular cloud models supports photo-excitation by the X-ray radi- 
ation from the active nucleus, which is not blocked by a molecular torus in this direction. 
However, with current data, the contribution from UV fluorescence or collisional excitation 
cannot be ruled out. 



^The X-rays from the AGN also explain the recent non-detection of CO gas (IDumas et al. 

2010 ) in contrast to the presence of 10^— lO^M© molecular hydrogen gas (IStorchi-Bergmann et al, 
2OO9I ) . The total H2 mass derived from the X-ray observation agrees with that expected from 
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the mass of "hot" H2. The absorption column derived from the X-ray spectral fit, when com- 
pared to the XDR models, implies that the lack of CO close to the nucleus may be explained 
by its low abundance relative to H2 in the XDR. 
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Fig. 1. — ACIS raw image of the NGC 4151 nuclear region (a) in the soft band (0.3-1 keV). 
(b) the hard band (1-7 keV). X-ray enhancements "A" and "B" are indicated with arrows. 
Sohd hnes mark the directions of the kinematic major axis of the host galaxy (P.A.~22°, 
Pedlar et al. 1992; Mundell et al. 1999), the large scale "weak fat bar" (P.A.~130°; Mundell 
& Shone 1999), and the ENLR bicone (P.A.~65°, Evans et al. 1993). The narrow streak 
running north-south across the nucleus (present in both panels) is an artifact of the CCD 
readout of the bright nucleus. 
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Fig. 2. — (a) Comparison of model PSFs (dashed line: without pile-up effect; dotted line: 
pile-up included) with the observed radial profile of the central bright point-like source in 
the hard band (1-7 keV). (b) Comparison of model PSFs with the observed radial profile in 
the 0.3-1 keV band for the obviously extended NE-SW sector, (c) the same as (b) but for 
the less clearly extended SE-NW sector, (d) Radial variation for the ratio between events 
with bad grades (1,5,7) and those with good grades (0,2,3,4,6). The level of such ratio in 
the Chandra ACIS-stowed background (pile- up free) is shown (dotted hne). Both the 0.3-1 
keV and 1-7 keV bands are shown. 
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(a) Smoothed ACIS 1/8 pixel 



(b) Smoothed HRC 
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Fig. 3.— (a) 0.3-1 keV ACIS subpixel image of NGC 4151 (1/8 native pixel); (b) 0.1- 
10 keV HRC image. Both are smoothed by a FWHM = 0.4" Gaussian kernel for better 
visualization of faint extended features. 
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(d) ACIS 0.3-1 keV (subpixel) 

Obsid 372 (0.1s frame) 
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Fig. 4.— (a) Raw 0.3-1 keV ACIS image of the inner 6" x 6" NGC 4151 nucleus; (b) 
ACIS image of the same region with subpixel binning (1/8 native pixel), demonstrating the 
improved resolution; (c) Same as (b) but using only single pixel events (grade 0); (d) 0.3-1 
keV ACIS image (0.1 s frame time) from ObsID 372 with subpixel binning; (e) same as (d) 
but for a frame time of 0.4 s; (f) HRC image (0.1-10 keV). All subpixel resolution image were 
smoothed with a FWHM = 0.25" Gaussian kernel for better visualization of faint extended 
features. 
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Fig. 5.— (a) PSF-deconvolved ACIS image (0.3-1 keV) using the EMC2 method (100 
iterations; Esch et ah 2004; Karovska et al. 2007); (b) PSF-deconvolved HRC image (0.1-10 
keV; Wang et al. 2009) using the EMC2 method. Contours of the deconvolved HRC image 
shown in (b) are overlaid in both panels to assist comparison of the features. 



- 28 - 




0.002 0.0059 0.014 0.03 



Fig. 6. — (a) Hardness ratio map of the central 1' x 1' (3.9 kpc-across) region, defined as the 
ratio of the counts in the 1-2 keV to the 0.3-1.0 keV band. Contours in white (sohd) and 
red (dashed) show the continuum-subtracted Ha (iKnapen et aLl I2004J ) and ^^CO emission 
(IDumas et al.ll2010l ). respectively, (b) Same as (a) but zoomed-in to the inner 6" region 
outlined by a box in (a). Th e cross marks the location of t he nucleus. Grey contours show 
the H2 A2.1218 /xm emission (jStorchi-Bergmann et al.ll2009l ). 
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Fig. 7.— V-H color map of the NGC 4151 nuclear region from HST F550M and F160W 
band images. Point source has been removed using the simulated PSF and the resolution of 
F550M image was degraded to match that of the F160W image. Contours (red) are from the 
X-ray hardness ratio map shown in Figure [6]d. Note that residuals from diffraction spikes 
can still be seen at P.A.^407220°. 




Fig. 8. — The spectral extraction regions overlaid on the hardness ratio image. The poly- 
gons outline the regions Rl (inner cone region), R2 (edges of the cone), and R3 (region 
perpendicular to the bicone-axis). 
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Fig. 9. — The X-ray spectra and fits obtained for three regions defined in Figure |9] that are 
morphologically distinct (see text and Table [2]). 
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Table 1. Log of Chandra ACIS-S Observations of NGC 4151 



Obs ID 


Dale- 


liisl luuinit 


IVaiiic TiuK! (s) 


Exi>cjsui(- (ks) 


372 


Mar 06, 2000 


ACIS-S 


0.1 


0.8 


372 


Mar 06, 2000 


ACIS-S 


0.4 


6.2 


9217 


Mar 29, 2008 


ACIS-S 


0.6 


116 


9218 


Mar 27, 2008 


ACIS-S 


0.6 


63 


9219 


Mar 02, 2008 


HRC-I 


1.6 X 10-5 


49 



Table 2. Spectral Fits to Extended X-ray Emission Features 



Region^ 




Nh 
[xl022 cm-2] 


Temperature 
[keV] 


PSF fraction'^ 

[%] 


/ d.o.f. 


Rl 
R2 
R3 


'^•"•^-0.02 

-0.27 ±0.03 
-0.04 ±0.02 


0.02=* 
70+0-45 

6.5±1.5 


0.08«;°^ 


0.5 
0.7 
2.2 


135/94 
66/75 
261/231 



^The regions and spectra are presented in Figure |9] Rl-the central inner part of the 
cone; R2-the edges of the cone; and R3-the high hardness ratio region perpendicular to 
the bicone-axis. 

''Hardness ratio is defined as (Ci_2koV - C'o.3-ikcv)/(Ci_2kcV + C'o.S-lkov)- 
•^The fraction of nuclear emission contributed to the extended emission. 
''Absorption fixed at Galactic column as required by the fit. 



